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1.  Introduction 

Quantifying  the  energy  content  of  the  coupled  magneto¬ 
sphere-ionosphere-thermosphere  (MIT)  system  during  geomag¬ 
netic  storms  is  a  long  sought  goal  of  space  weather  research.  The 
Dessler-Parker-Sckopke  (DPS)  theorem  argues  that  the  total 
energy  of  ring-current  (£RC)  particles  is  directly  proportional  to 
the  Dst  index  (cf.  review  by  Carovillano  and  Siscoe,  1973).  The 
Assimilative  Mapping  of  Ionospheric  Electrodynamics  (AMIE) 
technique  (Richmond  and  Kamide,  1988)  combines  data  from 
ground  and  satellite-borne  sensors  to  reconstruct  electric  fields  (c) 
and  current  (j)  patterns  at  high  latitudes.  In  its  simplest 
representation  energy  into  the  ionosphere-thermosphere  (IT) 
takes  the  form  of  Joule  heat  =  Yp  c2,  where  Yip  represents  the 
Pedersen  conductance  of  the  ionosphere.  In  fact,  energy  coupling 
is  much  more  complicated  than  simple  Joule  heating  leads  to  the 
formation  of  stormtime  wind  systems  (Richmond  and  Thayer, 
2000)  that  drive  disturbance  dynamos  (Blanc  and  Richmond, 
1980).  Disturbance  dynamos  in  turn  alter  prevailing  electric  field 
patterns  throughout  the  MIT  system.  Many  aspects  of  this  high-to- 
low  latitude  electrodynamic  coupling  are  now  captured  in 
computer  simulations  (Maruyama  et  al..  2005). 

Energy  enters  the  stormtime  ionosphere  at  auroral  and  polar 
cap  latitudes  in  the  form  of  particle  precipitation  and  Poynting 
flux  (Knipp  et  al..  2005)  and  is  mostly  transferred  to  the 
thermosphere  via  ion-neutral  collisions.  Consequent  heating  of 
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neutrals  alters  thermospheric  scale  heights  allowing  molecular 
species  access  to  higher  than  quiet-time  altitudes  and  generates 
winds  that  transport  energy  to  subauroral  latitudes.  Quantitative 
understanding  of  how  energy  transport  via  neutral  winds  still  has 
significant  gaps.  Bruinsma  and  Forbes  (2007)  identified  a 
sequence  of  traveling  atmospheric  disturbances  (TADs)  in  mid- 
to-low  latitude  measurements  by  the  accelerometer  on  the 
Challenging  Minisatellite  Payload  (CHAMP)  satellite  during  a 
magnetic  storm  with  multiple  heat  injections.  Molecular  scale 
heights,  altered  during  TAD  encounters,  produce  effects  that  are 
identifiable  through  remote-sensing  techniques.  Lin  et  al.  (2005) 
showed  that  during  the  October-November  2003  superstorm 
both  the  [0)/|N2l  ratios  and  total  electron  content  (TEC)  of  the 
equatorial  ionization  anomaly  decrease.  The  latter  reflects 
enhanced  recombination  rates  due  to  increased  molecular  con¬ 
centrations  in  the  F  layer. 

This  report  extends  our  preliminary  study  of  IT  responses 
observed  by  accelerometers  on  the  Gravity  Recovery  and  Climate 
Experiment  (GRACE)  satellites  during  the  magnetic  storms  of 
November  2004.  Burke  et  al.  (2007a)  found  that: 


(1)  Presently  used  models  under-predict  stormtime  density 
increases  by  >  100%, 

(2)  Local  densities  p  varied  widely  but  global  averages  f>  evolved 
systematically. 

(3)  Satellites  crossing  the  polar  cap  near  the  noon-midnight 
meridian  during  magnetic  storms  experience  positive  and 
negative  “density”  spikes.  These  spikes  are  head-  and  tail¬ 
wind  effects  due  to  coupling  with  anti-sunward  convecting 
ions. 
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(4)  Magnetosphenc  electric  fields  £vs  -  ^pc/2/.vRe  predicted  from 
parameters  observed  at  L,  anticipate  variations  of  p  with  lead 
times  of  4-5  h.  <Z>pc  and  2LyRe  are  the  polar  cap  potential  and 
width  of  the  magnetosphere,  respectively. 

(5)  During  the  main  and  early  recovery  phases  p  correlates  with 
Dst. 

(6)  The  thermosphere  relaxes  to  pre-storm  values  much  faster 
than  Dst. 

Allowing  for  local  turbulence,  result  (2)  suggests  that  on  a 
global  scale,  the  stormtime  thermosphere  evolves  as  a  large 
thermodynamic  system  that  never  strays  far  from  diffusive 
equilibrium  (Wilson  et  al.,  2006;  Burke,  2008).  The  present 
analysis  suggests  that  the  stormtime  thermosphere  responds  to 
interplanetary  forcing  as  a  driven-dissipative  system  whose 
dynamics  are  described  by  a  simple  differential  equation. 

In  this  paper,  we  compare  global  estimates  of  IT  energy 
budgets  during  the  magnetic  storm  of  November  2004.  Calcula¬ 
tions  employ  three  independent  approaches.  The  first,  involves 
calculating  net  (incident-reflected)  Poynting  flux  (S»)  into  the 
ionosphere  using  electric  (<5k)  and  magnetic  (<5B)  field  perturba¬ 
tions  measured  along  trajectories  of  Defense  Meteorological 
Satellite  Program  (DMSP)  spacecraft  (Huang  and  Burke,  2004). 
The  second,  a  statistical  model  developed  by  Weimer  (2005)  (W5), 
predicts  distributions  of  S  into  the  northern  and  southern  high- 
latitude  ionospheres  using  the  solar  wind  (SW)  and  interplanetary 
magnetic  field  (IMF)  as  drivers.  The  third  method  uses  orbit- 
averaged  densities  (p)  inferred  from  accelerometer  measurements 
on  the  polar-orbiting  GRACE  satellites  to  estimate  thermospheric 
energy  (Eth)  content  at  altitudes  above  100  km  (Burke,  2008).  The 
following  sections:  (1)  describe  the  methodologies  and  data 
sources,  (2)  provide  examples  of  energy  and  power  budgets 
during  the  November  2004  magnetic  storm  via  the  three  different 
methods,  and  (3)  comment  on  some  implications  of  resultant 
similarities  and  differences 

2.  Methodologies  and  data  sources 

The  first  method  estimates  S„  from  drift  meter  and  magnet¬ 
ometer  measurements  acquired  along  the  trajectories  of  DMSP 
spacecraft  as  they  cross  the  high-latitude  ionosphere  at  840  km. 
Using  a  transmission  line  analogy,  Huang  and  Burke  (2004) 
demonstrated  that  S  *{Sey  x  <5B7)/p0  is  the  net  (incident-re¬ 
flected)  Poynting  flux.  In  making  such  calculations  no  knowledge 
about  the  distribution  of  conductivity  beneath  the  spacecraft  is 
required.  Baselines  for  electric  and  magnetic-field  calculations 
are  established  using  the  offset-subtraction  technique  described 
by  Rich  et  al.  (2007).  Because  it  would  be  necessary  to 
approximate  energy  depositions  far  from  spacecraft  trajectories, 
this  is  not  a  practical  tool  for  making  realistic  estimates  of  total 
power/energy  inputs  to  the  global  IT  system.  However,  having 
direct  measures  of  the  actual  electromagnetic  energy  inputs  along 
trajectory  lines  provides  quantitative  means  for  evaluating  and 
perhaps  improving  the  predictions  of  global  models  such  as  W5  or 
AMIE. 

Weimer  (2005)  combined  two  statistical  models,  initially 
developed  to  estimate  the  distributions  of  electric  potential  <P{A, 
<l>)  (Weimer,  1995)  and  field-aligned  currents  j  (/l,  (f> )  (Weimer, 
2001)  in  the  high-latitude  ionosphere,  to  determine  S„  as  a 
function  of  magnetic  latitude  /l,  magnetic  local  time  (p  as  well  as 
the  interplanetary  magnetic-field  magnitude  and  orientation. 
Maps  of  Poynting  flux  (S  ),  produced  at  5  min  intervals,  reflect 
evolving  solar-wind  and  IMF  conditions.  Total  input  power  (PWs) 
is  determined  by  integrating  S  over  affected  areas  of  the 
ionosphere.  Applying  the  first  law  of  thermodynamics,  Wilson 


et  al.  (2006)  tested  W5  by  comparing  quiet-time  mass  densities 
inferred  from  the  drag  exerted  on  75  reference  satellites  at 
altitudes  between  200  and  800  km  (Casali  et  al.,  2002;  Storz  et  al., 
2002)  with  their  disturbed-time  counterparts.  Increased 
densities  observed  during  periods  of  high  magnetic  activity  were 
attributed  to  increases  in  the  thermal  energy  content  of  the 
thermosphere  and  work  done  against  the  force  of  gravity.  Energy 
increases  inferred  from  the  reference  satellite  data  were  similar  to 
but  less  than  predictions  of  W5,  Since  some  of  fraction  of  the 
energy  must  be  deposited  at  altitudes  below  200  km,  these 
differences  were  regarded  as  being  consistent  with  the  model’s 
general  correctness. 

The  third  approach  uses  orbit-averaged  mass  densities  p 
measured  by  accelerometers  on  the  polar-orbiting  GRACE  satel¬ 
lites  (Tapley  et  al.,  2004)  to  infer  the  total  energy  £th  of  the 
thermosphere  at  altitudes  /i>100km.  Because  the  technique  is 
relatively  new,  a  fuller  explanation  is  required.  As  indicated  in 
Fig.  1  the  Jacchia  model  (Jacchia  1977)  (J77)  consists  of  three 
elements.  The  first  is  a  set  of  empirically  based  driving  equations 
that  relate  variation  of  solar  fluxes  and  the  ap  index  to 
distributions  of  exospheric  temperatures  T ,  (Bowman  et  al., 
2006).  The  second  consists  of  15  tables  that  describe  the 
thermospheric  density  and  temperature  T  profiles  required  for 
diffusive  equilibrium  with  given  values  of  T The  tables  list 
profiles  for  700 ^ 2000 K  in  100 K  increments.  The  third 
element  consists  of  predicted  densities  and  consequent 
atmospheric  drag  exerted  on  space  objects. 

Burke  (2008)  argued  that  since  T*  uniquely  specifies  each 
p-T-h  profile  in  Jacchia  tables  and  GRACE  measures  p  at  known 
altitudes,  it  should  be  possible  to  drive  the  model  directly  by 
accelerometer  measurements.  Treating  table  elements  of  the  J77 
model  as  though  they  were  ’Teal  data,"  statistical  analyses 
revealed  that  to  a  good  approximation 

2 

T-»  =  ]Ta((/i)p'(/i>  (1) 

i-0 


Fig.  1.  Flow  chart  contrasting  standard  procedures  (left)  for  estimating  the 
atmospheric  drag  exerted  on  objects  flying  in  low  Earth  orbit  from  observed  solar 
UV  radiances,  and  the  ap  index  with  a  new  method  (right)  using  orbit-averaged 
mass  densities  inferred  from  measurements  of  accelerometers  an  the  polar- 
orbiting  CHAMP  and  GRACE  satellites. 
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Over  the  range  300<h<500km  each  coefficient  a,<h)  is  well  SS1ES  consists  of  spherical  Langmuir  probes  on  0,8  m  booms  to 

represented  by  a  fifth-order  polynomial  in  h.  measure  electron  densities  and  temperatures  and  three  different 
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Regression  coefficients  obtained  in  fitting  Tx  and  a*  to 
polynomials  exceed  0.999.  In  applying  Eqs.  (1)  and  (2)  /)  and  h 
are  in  grams/cc  and  kilometers,  respectively.  During  the  period  of 
interest  the  orbit-averaged  altitude  of  GRACE  was  -486.7  km. 

Using  information  contained  within  J77  tables  and  following 
the  method  of  Wilson  et  al.  (2006),  Burke  (2008)  demonstrated 
that  the  total  (thermal  and  gravitational)  energy  of  the  thermo¬ 
sphere,  in  Joules,  at  h>100km  is  linearly  related  to  the  globally 
averaged  exospheric  temperature  ty 

=  5.365  x  1017  +  8  727  x  10I3T„  (3) 

Over  the  Tx  range  700-2000  K,  Eth  only  changes  by  -15%. 
Experience  gained  analyzing  three  years  of  data  from  the  CHAMP 
and  GRACE  satellites  indicates  that  during  most  storms  Tr 
increased  by  a  few  hundred  K.  In  comparison  to  base  lines  that 
reflect  prevailing  solar  ultraviolet  (UV)  driving,  stormtime 
percentage  increases  of  Eth  are  small. 

Burke  (2008)  noted  two  other  relevant  properties  of  the  />,  T  . 
Eth  triad:  (1)  they  vary  on  time  scales  of  solar  rotations  and 
geomagnetic  activity  represented  by  and  the  Dst  index.  (2) 
Because  they  are  related  linearly,  Eth  and  t ,  decrease  at  the  same 
rates  whenever  sVs  turned  off.  Plots  of  lnEth  versus  time  (not 
shown)  during  magnetic  storms  indicate  that  as  Cvs  —  0.  E,h  relaxes 
with  an  e-fold  time  T^6.5h.  The  first  property  allows  us  to  treat 
Eth  as  derived  from  two  independent  but  additive  sources,  solar 
ultraviolet  radiation  Eth  uv  and  the  solar  wind  Eth  Sw-  As  indicated 
below,  the  second  property  allows  us  to  use  time  histories  of  r,v s  to 
estimate  Eth  sw  and  then  to  compare  its  temporal  variations  with 
those  determined  from  GRACE  measurements  of  p  with  power 
input  predictions  of  W5. 

2.1  Data  sources 

The  solar-wind  density  and  velocity  were  measured  by  the 
solar  wind  electron,  proton,  and  alpha  monitor  (SWEPAM) 
(McComas  et  al.,  1998)  and  the  IMF  vector  was  observed  by  the 
magnetic-field  instrument  (MF1)  (Smith  et  al.,  1998)  on  the 
Advanced  Composition  Explorer  (ACE)  spacecraft  in  a  halo  orbit 
around  L\,  -234  RF  upstream  of  Earth. 

Electric  and  magnetic-field  perturbations  were  measured  by 
ion  drift  meters  (1DM)  and  triaxial  fluxgate  magnetometers  on 
DMSP  FI 3.  F15,  and  F16.  DMSP  spacecrafts  are  three-axis 
stabilized  and  fly  in  circular,  sun-synchronous,  polar  (inclination 
98.7  )  orbits  at  an  altitude  of  -840  km.  The  geographic  local  times 
of  the  orbital  planes  are  near  the  18:00-06:00  (F13),  20  00-08:00 
(F16).  and  2100-0900  (FI 5)  LT  meridians.  Offsets  between  the 
geographic  and  geomagnetic  poles  allow  DMSP  spacecraft  to 
sample  wide  ranges  of  magnetic  local  times  over  the  course  of  a 
day.  Each  satellite  carries  a  suite  of  sensors  that  measures  the 
densities,  temperatures,  and  drift  motions  of  ionospheric  ions  and 
electrons  (Special  Sensor  for  Ions  Electrons  and  Scintillations, 
SS1ES),  and  perturbations  of  the  Earth’s  magnetic  field  (Special 
Sensor  for  Magnetic  Fields.  SSM). 


sensors  mounted  on  a  conducting  plate  facing  the  satellite  ram 
direction.  They  are:  (1)  an  ion  trap  to  measure  the  total  ion 
densities;  (2)  an  1DM  to  measure  horizontal  (VH)  and  vertical  (Vv) 
cross-track  components  of  plasma  drifts;  (3)  a  retarding  potential 
analyzer  (RPA)  to  measure  ion  temperatures  and  the  in-track 
component  of  plasma  drift  Vn  (Rich  and  Hairston,  1994)  Relatively 
small  uncertainties  about  prevailing  spacecraft  potentials  lead  to 
unacceptable  errors  in  V  calculations  and  are  not  included  in  this 
study.  Thus.  Poynting  fluxes  calculated  from  DMSP  measurements 
are  lower-bound  estimates.  VH  and  Vv  are  sampled  6  times  per 
second  but  are  normally  reported  as  4  s  averages.  They  are  then 
used  to  calculate  the  in-track  component  of  electric  field  (£y0  and 
the  distribution  of  electric  potential  (</>)  along  the  trajectory. 
DMSP  estimates  of  the  polar  cap  potential  (<J»pc)  are  reported  as 
differences  between  measured  potential  extrema  in  the  northern 
and  southern  hemispheres.  Since  trajectories  usually  do  not  cross 
potential  maxima  and  minima,  we  also  regard  DMSP  as  providing 
lower-bound  estimates  of  the  true  </>po 

The  triaxial  fluxgate  SSM  sensors  are  either  mounted  on  the 
spacecraft  body  (F13)  or  on  5m  booms  (F15,  F16)  to  reduce 
susceptibility  to  spacecraft-generated  contamination.  Magnetic 
vectors  are  sampled  16s“l  and  reported  as  Is  averages  of 
SB  =  Bmeas-BiC(Rh  differences  between  measured  and  Interna¬ 
tional  Geomagnetic  Reference  Field  (1GRF)  values  of  magnetic 
fields  at  the  spacecraft  locations.  Data  are  presented  as  SBX,  SBY, 
SB 7  in  a  spacecraft-centered  coordinate  system  TheX-  and  V-axis 
point  toward  spacecraft  nadir  and  along  its  velocity  vector, 
respectively.  The  Z-axis  completes  the  right-hand  system  and 
generally  points  in  the  anti-sunward  direction.  The  downward 
component  of  Poynting  flux  is  estimated  as  S  *(&:y  x  SBz)frto. 
Values  of  S  presented  below  were  calculated  and  reported  at  a 
cadence  of  one  sample  per  second. 

Neutral  densities  at  —486.7  km  were  estimated  from  measure¬ 
ments  of  Spatial  Triaxial  Accelerometer  for  Research  (STAR) 
sensors  on  two  GRACE  satellites  flying  in  tandem  in  nearly 
circular,  polar  orbits  (Tapley  et  al.,  2004).  STAR  sensors  monitor 
electrostatic  forces  needed  to  maintain  proof  masses  (PM)  at  the 
center  of  a  cage  located  within  2  mm  of  each  spacecraft’s  center  of 
mass  (Bruinsma  et  al.,  2004).  Both  the  spacecraft  and  PM  respond 
to  gravity  in  the  same  way.  Thus,  changes  in  electrostatic  forces 
that  maintain  the  PM  at  its  cage’s  center  reflect  spacecraft 
responses  to  non-gravitational  forces  such  as  atmospheric  drag 
and  radiation  pressure  (Bruinsma  and  Biancale,  2003).  Accelera¬ 
tion  due  to  drag  is  given  by 

fldr.j  =  CdA!C/2Mx(>V2  (4) 

where  Asc  and  Msc  represent  the  effective  cross-sectional  area  and 
mass  of  the  spacecraft,  respectively,  p  the  mass  density  of  the 
neutral  atmosphere,  and  V  the  spacecraft  velocity  in  the  rest  frame 
of  ambient  neutrals.  The  drag  coefficient  C0  depends  on  the  angle 
of  flow  to  the  spacecraft  surface,  the  ratio  of  the  temperatures  of 
the  satellite  surface  and  the  local  atmosphere,  and  the  ratio  of  the 


312 


W.J.  Burke  et  aL  /  Journal  of  Atmospheric  and  Solar-Terrestrial  Physics  72  (2010)  309-318 


mean  mass  of  atoms  in  the  atmosphere  to  those  on  the  satellite 
surface  (Bruinsma  and  Biancale,  2003). 


3.  Observations 

This  report  focuses  on  electrodynamic  responses  observed  on 
November  10,  2004.  Here,  we  refer  to  both  calendar  dates  and 
Julian  days  (JD)  of  the  year,  where  January  1  corresponds  toJD  1 
Reported  events  of  JD  315  occurred  within  a  larger  disturbance. 
Fig.  2  provides  context  in  the  form  of  the  Dst  index  (top),  the  ring- 
current  energy  ERC  (middle),  and  £Vs  (bottom)  for  JD  310-316 
(November  5-11).  ERC  was  calculated  using  a  form  of  the  DPS 
theorem  given  by  Stem  (2005):  £rc-(3  Dst  £M/2B0),  where 
Bo-3.1xl04nT  is  the  magnetic  field  at  the  equator,  and 
£M5:8x10,7J  the  magnetic-field  energy  above  the  Earth’s 
surface.  Methods  used  to  calculate  £Vs  are  described  in 
Appendix  A. 


Plots  in  Fig.  2  show  that  the  disturbance  began  late  on  JD  312 
(November  7)  after  a  long  quiet  period.  The  first  main  phase 
reached  a  maximum  epoch  at  06:00  UT  on  JD  313  (November  8) 
when  Dst  and  ERC  reached  -372  nT  and  -1.5  x  1016J,  respectively. 
Figures  3d  and  4d  of  Burke  et  al.  (2007b)  show  that  the  quantity 
/vs  =  /cvs(t)dt  highly  correlates  with  Dst  during  the  main-phase 
episode  of  JD  313,  Magnetic  quieting  was  interrupted  late  on  JD 
314  by  a  brief,  but  intense,  southward  turning  of  the  IMF  that  led 
to  a  spike  in  tVs  of  ^2.3  mV/m  near  19:00  UT.  The  geomagnetic 
disturbance  resumed  on  JD  315  with  Dst  reaching  a  minimum  of 
-289  nT  (£rc*  1.2  x  1016J)  at  09:00  UT. 

Plots  in  Fig.  3  show  thermospheric  responses  as  observed  by 
GRACE  in  polar  orbit  near  the  noon-midnight  meridian.  The  figure 
shows  orbit-averaged  densities  (top),  exospheric  temperatures 
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Fig.  2.  History  of  the  Dst  index  (top)  from  November  S— 11.  2004  (Julian  days 
310-316),  the  energy  of  the  ring-current  £RC  estimated  via  the  Dessler-Parker-Sck- 
opke  relation  (middle)  and  the  Volland-Stern  electric  field  cv s  determined  from 
ACE  measurements  near  Lt  (bottom)  plotted  as  functions  of  UT. 


Fig.  3.  GRACE  measurements  of  orbit-averaged  mass  densities  (top),  exospheric 
temperatures  (second),  and  thermospheric  energy  (third)  plotted  as  functions  of 
UT  for  (November  6-14)  310-317  (November  6-14).  2004.  The  Kvs  trace  is  repeated 
for  reference  in  the  bottom  panel. 
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(second),  total  thermospheric  energies,  and  eVs  (bottom)  plotted 
as  functions  of  time  over  the  same  7-day  interval  as  in  Fig.  2.  T , 
and  Eth  were  calculated  using  Eqs.  (1)  and  (3).  respectively. 
Attention  is  directed  to  six  empirical  features  apparent  in  the 
data: 

(1)  Decreasing  trends  in  traces  of  p,  fx,  and  Eth  on  JD  310  and  311 
when  f>vs  was  small  are  solar  rotation  effects  indicated  by  a 
decreasing  F}07  (not  shown). 

(2)  At  the  beginning  of  JD  312,  t ,  and  Eth  were  ~925K  and 
6.15  x  017J,  respectively. 

(3)  t r.  and  Eth  maxima  reached  ~1400K  and  6.6  x  1017J. 
respectively  on  JD  313. 

(4)  E{h  increased  by  5.5  x  1016,  while  ERt  only  reached  1.2  x  1016 
on  JD  315. 

(5)  T,  and  Eth  attained  maxima  of  ~1350K  and  6.55  x  1017 J, 
respectively  on  JD  315. 

(6)  T ,  and  Et h  fell  at  the  same  rate  on  JD  313  and  315  after  cVs 
turned  off. 

Fig.  4  illustrates  the  procedure  used  for  comparing  W5 
predictions  of  S  with  DMSP  measurements.  The  left  plot  shows 
a  typical  W5  Poynting  flux  map  on  a  (/1-0)  grid.  Iso-contour  lines 
indicate  the  spatial  distribution  of  electromagnetic  power  into  the 
southern  hemisphere  near  the  end  of  the  main  phase  (09:45  UT) 
of  the  November  10  disturbance.  The  W5  map  shows  two  local 
maxima  in  S  at  stormtime  auroral  latitudes  of  58  and  62  MLat 
near  dusk  and  dawn,  respectively.  Significant  but  lesser  power 
flows  into  the  polar  ionosphere  Curiously,  the  model  assigns 
almost  no  electromagnetic  power  input  to  the  regions  of  the 
dayside  cusp  or  the  midnight  sector  of  the  auroral  oval.  Indicated 
to  the  lower  left  and  right  of  the  map  are  magnitudes  of  minimum 
and  maximum  S .  Integrated  power  into  the  southern  polar  cap 
(886  GW)  is  provided  in  the  upper  right  corner.  The  solid  orange 
line  across  the  contour  map  indicates  the  dawn-to-dusk  trajectory 
of  DMSP  FI 3  across  the  southern  polar  cap  at  about  the  same  UT. 
The  trajectory  passed  through  the  predicted  morning  side 
maximum,  but  missed  the  one  on  the  evening  side  The  plot  on 
the  right  side  of  Fig.  4  compares  W5  Poynting  flux  predictions 
along  the  F13  orbit  with  data  from  combined  1DM  and  SSM 
measurements.  With  obvious  caveats,  this  instance  shows 
reasonable  agreement  between  W5  predictions  and  DMSP 
observations.  However,  most  W5  and  DMSP  S  distributions  do 
not  agree  so  closely.  We  note  that  as  FI 5  crossed  the  southern 
auroral  oval  at  ~  10:03  UT  it  recorded  peak  SM  values  of  -40mW/m2 


near  MLat  =  -57  ,  MLT  =  23.0,  where  the  W5  map  assigns  very 
little  Poynting  flux. 

To  support  quantitative  comparisons  between  the  thermo¬ 
spheric  power  inputs  predicted  by  W5  and  measured  by  DMSP, 
without  becoming  bogged  down  with  local  variations,  we  decided 
to  compare  two  global  parameters.  W5  routinely  calculates  power 
into  the  northern  and  southern  ionospheres  every  5  min.  We  add 
these  together  to  obtain  the  predicted  global  Pws.  We  also 
integrated  S  (W/m)  measurements  along  full  trajectories  of 
DMSP  F15  and  F16  to  calculate  a  quantity 

n3(W/m)  =  <J>  S|  ds  =  Vjjt  /  S||  dt  **  7.5  x  1  O'*  ^  (OEy  x  <)Bz)t/ jig 

(5) 

The  absence  of  {6EZ  x  dBY\ljig  from  S„  in  Eq.  (5)  reflects  our 
inability  to  extract  reliable  estimates  of  the  in-track  component  of 
plasma  drift  from  current-voltage  sweeps  of  retarding  potential 
analyzers  on  DMSP  satellites  (Rich  and  Hairston,  1994).  Integrals/ 
sums  extend  from  one  ascending  or  descending  node  to  the  next 
and  are  reported  at  the  times  of  the  second  equatorial  crossing. 
The  top  panel  of  Fig.  5  contains  results  of  S  calculations  for  DMSP 
F15  (diamonds)  and  F16  (circles)  on  November  9  and  10.  The  plot 
in  the  middle  shows  the  history  of  global  power  input  predicted 
by  W5.  For  reference  £Vs  is  shown  at  the  bottom  Data  presented  in 
Fig.  5  support  three  empirical  conclusions: 

(1)  Agreement  between  FI  5  and  F16  measurements  of  S  was 
close  on  JD  314.  While  agreement  remained  good  on  JD  315, 
the  correlation  clearly  degenerated.  This  is  primarily  a  spatial 
effect.  On  both  days  the  separation  of  ascending  nodes  for  the 
two  spacecraft  was  —3  min  in  UT  and  ~-l  h  in  LT.  A  linear 
correlation  analysis  (not  shown)  indicates  that  [J 
F1S  =  1.85+0.91  []  F)6;  the  correlation  coefficient  is  ~0.83. 
The  <  1  slope  agrees  with  the  sense  of  W5  zonal  gradients 

(2)  While  it  is  possible  to  establish  correspondence  between  [\ 
structures  and  events  in  the  Pw s  trace,  overall  agreement  does 
not  appear  strong. 

(3)  However,  remarkable  agreement  is  evident  between  Pws  and 
£vs  features. 

We  comment  on  the  significance  of  the  2nd  and  3rd  observa¬ 
tions  in  the  discussion  section. 

While  our  emphasis  has  been  on  global-scale  variations  of 
thermospheric  parameters,  we  recognize  that  this  is  not  the  whole 
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Fig.  4.  Contours  of  Poynting  flux  Sj|  predicted  by  WS  for  the  southern  hemisphere  at  09  4S  UT  on  November  10,  2004  (left).  The  dawn-to-dusk  trajectory  of  DMSP  FI 3  is 
indicated  by  an  orange  line  across  the  map.  Plolson  the  righl  side  compare  F13  Poynting  flux  measurements  (black)  with  WS  predictions  (red)  along  the  F13  trajectory.  (For 
interpretation  of  the  references  to  colour  in  Ihis  figure  legend,  ihe  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  5.  Orbil-integrated  Poynting  flux  S  (top)  measured  by  DMSP  F15  (diamonds) 
and  F16  (circles)  and  global  power  PWs  input  predicted  by  W5  (middle),  plotted  as 
functions  of  UT  for  JD  3 14  and  315  (November  10-11).  2004.  The  cv s  trace  (botlom) 
is  provided  for  comparison  with  Pw5. 


story.  Data  shown  in  Fig.  6  illustrate  some  of  the  complex 
dynamics  that  dominate  local  measurements.  The  traces  show 
neutral  densities  measured  by  GRACE  during  four  successive 
nighttime  orbital  passes  plotted  as  functions  of  geographic 
latitude.  In  November  2004,  the  satellite  trajectory  passed  from 
north  to  south  close  to  the  midnight  meridian.  Thus.  GRACE 
crossed  the  southern  polar  cap  from  night  to  day.  Above  the 
panels  we  indicate  the  universal  times  of  equatorial  crossings  on 
November  10  During  the  first  orbit  GRACE  detected  a  density 
maximum  slightly  north  of  the  equator.  Observed  densities  varied 
by  a  factor  of  ~2  across  all  nightside  latitudes.  In  the  second  panel, 
labeled  02:50:02  UT.  a  large  spike  appears  at  high  southern 
latitudes,  that  is  in  part  a  headwind  effect  (Burke  et  al.,  2007a), 
while  a  broad  featureless  valley  extended  across  the  middle  and 
low-latitude  thermosphere  By  the  time  of  the  equatorial  crossing 
at  04:24:27  UT,  three  distinct  density  maxima  were  seen  at 


latitudes  extending  from  -80  to  60  .  By  the  time  of  the  05:58:17 
UT  pass,  the  highest  densities  were  once  more  found  near  the 
equator.  Densities  measured  in  these  maxima  are  about  three 
times  larger  than  those  observed  at  mid-to-high  latitudes.  We 
discuss  possible  connections  between  these  observations  and 
TADs  reported  by  Bruinsma  and  ^orbes  (2007)  in  the  following 
section. 


4.  Discussion 

In  the  previous  section,  we  applied  a  new  method  for 
estimating  the  thermosphere’s  total  energy  content  £,h  using 
orbital  averages  of  densities  inferred  from  measurements  of 
accelerometers  on  the  polar-orbiting  GRACE  satellites.  Data  were 
also  provided  that  allow  comparisons  (1)  between  the  estimated 
£lh  and  £RC  during  the  November  2004  storm,  and  (2)  with 
thermospheric  energy  and  power  requirements  inferred  from 
GRACE  measurements  and  W5,  respectively.  Based  on  these 
comparisons  and  simultaneous  DMSP  measurements  of  S„  we 
comment  on  the  Poynting  flux  distributions  of  the  W5  model. 
Finally,  we  briefly  discuss  the  mutually  supportive  roles  of  in  situ 
measurements  from  GRACE  and  information  gained  through 
remote  sensing  of  |0]/[N2]  ratios  and  the  TEC  content  of  the 
low-latitude  ionosphere. 

Huang  and  Burke  (2004)  reported  repeated  episodes  of  intense 
Poynting  flux  infusions  to  the  auroral  oval  during  the  main  phase 
of  large  magnetic  storms.  Four  DMSP  satellites  detected  these 
phenomena  over  a  wide  range  of  nightside  local  times.  This 
observed  spatial  distribution  suggested  that  energy  transferred  to 
the  magnetosphere  from  the  enhanced  solar  wind  largely 
accumulates  in  the  ring  current.  The  ring  current  would  thus 
function  as  a  reservoir  from  which  particle  and  electromagnetic 
energy  is  drawn  to  heat  the  stormtime  IT  system.  The  W5  map  in 
Fig.  4  appears  compatible  with  this  perception,  placing  most  to 
the  ionosphere  at  auroral  rather  than  polar  cap  latitudes. 

If  this  conjecture  is  correct,  we  expect  that  £RC  would  exceed 
stormtime  increases  in  £lh,  perhaps  by  large  amounts.  Calcula¬ 
tions  plotted  in  Figs.  2  and  3  indicate  the  opposite  is  true,  A£th^3 
£rc.  A  similar  relation  was  found  during  the  magnetic  storms  of 
July  2004  (Burke,  2008).  The  DPS  theorem  thus  indicates  that  the 
ring  current  cannot  act  as  a  reservoir  for  most  of  the  energy 
reaching  the  stormtime  thermosphere  as  suggested  by  Huang  and 
Burke  (2004).  Magnitudes  of  <Ppc  calculated  from  ACE  data  (cf. 
Appendix  A)  and  observed  by  DMSP  were  >200kV  on  both 
November  8  and  10.  2004.  To  supply  the  entire  peak  power  of 
~2-3TW  predicted  by  W5  (Fig.  5)  on  JD  315,  Pedersen  currents  of 
^5MA  would  have  to  flow  across  both  polar  caps  Since  some 
significant  fraction  of  the  total  power  must  enter  the  /-T at  auroral 
latitudes,  actual  trans-polar  current  demands  are  correspondingly 
less  than  the  5  MA  estimate. 

To  make  meaningful  comparisons  between  the  W5  predictions 
and  GRACE  measurements  it  is  critical  to  recognize  fundamental 
differences  between  them.  W5  uses  solar-wind/lMF  measure¬ 
ments  to  predict  distributions  of  Poynting  flux  into  the  northern 
and  southern  ionospheres.  Integration  over  the  affected  areas 
provides  the  global  rates  of  energy  input  On  the  other  hand,  orbit- 
averaged  £,h  measurements  from  GRACE  are  snapshots  that 
capture  the  total  energy  content  of  the  thermosphere  during  a 
given  interval  of  time.  They  reflect  thermospheric  give  and  take 
between  power  received  from  the  outside  and  lost  via  radiative 
and  other  processes  (Mlynczak  et  al.,  2005). 

We  noted  above  that  £th  varies  on  time  scales  that  correspond 
to  the  solar  rotation  period  and  changes  in  £vs  (Burke,  2008).  We 
have  found  that  as  £VS-*0.  £th  sw  relaxes  with  an  e-fold  time  scale 
of  r^6.5h.  This  behavior  suggests  that  the  global  thermosphere 
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Fig.  6.  Thermospheric  mass  densities  measured  by  GRACE  during  four  consecutive  passes  on  November  10.  2004  along  the  midnight  meridian  plotted  as  functions  of 
latitude.  The  universal  times  of  corresponding  equatorial  crossings  are  noted  above  the  plots. 
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Fig-  7.  Plots  of  t.vs  (black  line)  and  values  of  £th  sw  predicted  by  Eq.  (5)  (blue  line) 
and  measured  by  GRACE  (red  dots).  (For  interpretation  of  the  references  to  colour 
in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article  ) 


acts  as  a  driven-dissipative  system.  The  well-known  differential 
equation  describing  energy  changes  in  such  systems  has  the  form 

d^thsw  £thSW  ,a>. 

__=atvs___  (6) 

Similar  equations  govern  the  evolution  of  energy  in  R-L  circuits 
and  the  ring  current  during  magnetic  storms  (Burton  et  al.,  1975). 
Coupling  and  relaxation  constants  x  and  x  must  be  determined 
empirically.  Eq.  (6)  is  solved  numerically  by 

£ihswtfn+i)  =  ac.vsrtrjAt  4-  £thsw(Li)^l  — — ^  (7) 

Using  solar-wind  and  IMF  measurements  from  the  ACE  satellite 
we  calculated  hourly  averages  of  r.Vs-  For  convenience  we  take  the 
time  interval  At  to  be  1  h.  A  preliminary  analysis  of  GRACE  and  cvs 


Fig-  8.  Comparison  of  ou:vs  required  by  Eq,  (6)  to  match  GRACE  measurements  of 
with  PW5  inputs  predicted  for  November  9  and  10.  2004. 


measurements  indicates  that  if  we  consider  Eth  in  units  of  10lbJ, 
the  coupling  constant  a  ^0.55.  Fig.  7  compares  predictions  of 
Eq.  (7)  (blue  line)  with  values  of  £th  sw  inferred  from  GRACE 
accelerometer  measurements  (red  dots)  during  the  November 
storm  period.  The  black  line  repeats  the  history  of  the  *;Vs  driver 
Although  we  have  not  tried  to  optimize  estimates  of  a  and  r. 
agreement  between  predictions  and  measurements  appears  to  be 
close. 

Eq.  (6)  calculates  the  net  power  into  the  thermosphere  as 
power  input  (aeVs)  minus  the  dissipation  (£th  sw/t).  Thus,  xr.V5 
represents  the  power  input  needed  to  replicate  GRACE  measure¬ 
ments  and  should  be  directly  comparable  with  Cws.  If  we 
represent  £,„  in  Joules.  i«5.5x  10ls((J/h)/(mV/m))  =  1.528 TW/ 
(mV/m).  The  black  trace  in  Fig.  8  replicates  the  middle  panel  of 
Fig.  5,  showing  net  power  Pw s  TW;  red  dots  indicate 
corresponding  values  of  2£Vs  time-shifted  by  1  h  to  allow  for 
transport  from  L|  to  the  magnetosphere.  Correspondence  between 
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Pws  and  a£vs  appears  close.  Both  methods  indicate  that  at  several 
times  peak  power  inputs  reached  the  2-3 TW  range.  This 
significantly  exceeds  ~500CW  of  power  incident  to  the  dayside 
of  Earth  at  EUV  wavelengths  (Knipp  et  al.,  2005).  Agreement 
between  Pws  and  acv s  independently  confirms  the  J77  analysis  of 
GRACE  measurements  and  indicates  that  Pws  can  be  used  to  drive 
Eq.  (6)  to  predict  the  development  calculations  of  £th  sw  and  Trj. 

Siscoe  (1983)  demonstrated  that  only  “obliquely  propagating” 
Alfven  waves  carry  field-aligned  current.  These  waves  also  carry 
field-aligned  Poynting  flux  from  regions  where  V  5>0.  Both  J 
and  S  propagate  along  field  lines  that,  in  an  Earth  stationary 
frame  of  reference,  move  in  the  VD  =  (Ec  x  B)/B2  direction,  where 
Ec  represents  the  background  convection  electric  field.  Strictly 
speaking,  at  ionospheric  altitudes  the  field  line  is  not  exactly 
conjugate  to  the  J  and  S  source  locations.  However,  they  can  only 
be  removed  from  the  source  location  by  a  distance  VD-  AT*.  Here, 
ATa  is  the  Alfven  travel  time  required  for  the  signal  to  propagate 
from  the  magnetospheric  source  to  the  ionosphere.  At  typical 
plasma  sheet  and  ring-current  distances  the  separations  are 
relatively  small.  Data  in  Fig.  4  place  the  bulk  of  S  predicted  by  W5 
at  auroral  latitudes.  As  such  S  sources  appear  to  lie  in  the  plasma 
sheet  and  ring  current.  Data  presented  in  Figs.  2  and  3 
demonstrate  that  Eth  sw-3£RC.  Given  the  hour-long  time  scales 
needed  to  transport  and  energize  ring-current  particles,  it  is 
difficult  to  understand  how  energy  would  be  processed  in  the 
magnetosphere  in  ways  that  explain  the  observed  high  and 
sustained  rates  of  power  transported  to  the  stormtime  auroral 
ionosphere-thermosphere.  Resolution  to  this  apparent  discre¬ 
pancy  lies  in  the  fact  that  the  peak  magnitude  of  S*  is  located  at 
the  boundary  between  Regions  1  and  2  currents.  The  dividing  line 
between  the  solar-wind  and  ring-current  sources  of  Poynting  flux 
lies  somewhat  poleward  of  the  peak,  but  still  within  the  auroral 
oval  It  is  possible  that  an  area  integration  of  5  within  different 
regions  will  confirm  the  three  to  one  ratio,  but  this  calculation  is 
beyond  the  scope  of  this  observational  report. 

A  comparison  with  the  Dst  trace  in  Fig.  2  shows  that  latitudinal 
density  profiles  in  Fig.  6  were  acquired  as  magnetic  activity 
emerged  from  a  6  h  lull  (314.22  to  315.04  UT)  into  the  main-phase 
intensification.  The  intensities  of  peak  latitudinal  variations 
observed  near  04:24  and  05:58  UT  exceeded  those  of  troughs  by 
>50%.  We  tentatively  interpret  these  structures  as  signatures  of 
equatorward  propagating  TADs  similar  to  those  reported  by 
Bruinsma  and  Forbes  (2007). 

TADs  provide  evidence  of  energy  transport  from  the  auroral 
zones  where  heating  occurs  towards  the  equator  where  the 
atmospheric  pressure  is  lower.  However,  the  Earth  rotates  by  -23 
in  longitude  from  one  ascending-node  crossing  of  GRACE’S  orbit  to 
the  next.  To  interpret  the  structures  such  as  those  detected  by 
GRACE  (Fig.  6)  as  TADs,  it  is  necessary  to  assume  that  they  extend 
over  wide  swaths  of  longitude  (Bruinsma  and  Forbes,  2007).  Here, 
information  gained  through  remote-sensing  techniques  is  suppor¬ 
tive.  Because  of  their  effects  on  [0]/[N2]  ratios  and  TEC  at  low 
latitudes,  TAD  effects  can  be  associated  with  Composition  Perturba¬ 
tion  Zones  (CPZ)  (Lin  et  al.,  2005).  During  the  same  interval  as  the 
data  shown  in  Fig.  6.  Hei  (private  communication,  2008)  observed 
TEC  decreases  that  progressed  from  high  latitudes  in  the  southern 
(summer)  hemisphere  toward  the  equator.  Critical  for  the  inter¬ 
pretation  of  structures  detected  by  accelerometers,  Hei  and  co- 
workers  determined  from  UV  images  that  CPZs  do  indeed  maintain 
coherence  over  several  hours  in  local  time  (longitude). 


5.  Summary  and  conclusions 

In  this  paper,  we  employed  three  independent  techniques  to 
estimate  the  amount  of  energy  that  entered  the  high-latitude 


ionosphere-thermosphere  as  Poynting  flux  during  the  magnetic 
storm  of  November  2004.  (1)  Data  from  plasma  drift  and 
magnetic-field  sensors  on  the  DMSP  F15  and  F16  spacecraft 
combine  to  support  calculations  of  Poynting  flux  along  trajec¬ 
tories.  (2)  The  empirical  model  of  Weimer  (2005)  uses  plasma  and 
magnetic-field  measurements  acquired  near  to  predict  Poynt¬ 
ing  flux  distributions  in  magnetic  latitude  and  magnetic  local  time 
at  a  5-min  cadence.  Total  power  Pws  input  is  obtained  by 
integrating  over  the  affected  area  of  the  upper  atmosphere.  (3) 
We  introduced  a  new  technique  that  uses  orbit-averaged  mass 
densities  measured  by  the  polar-orbiting  GRACE  satellites  to  infer 
Eth  the  total  energy  content  of  the  thermosphere.  Values  of  £th 
above  quiet-time  levels  were  compared  with  £RC  the  energy 
content  of  the  ring  current  estimated  via  the  Dessler-Parker-Sck- 
opke  relation.  Observed  properties  of  Eth  responses  to  eVs  driving 
indicate  that  the  global  thermosphere  acts  as  a  driven-dissipative 
system  whose  behavior  is  described  by  a  simple  first-order 
differential  equation.  The  coupling  (a)  and  relaxation  (t)  constants 
of  the  equation  were  established  empirically  through  comparisons 
with  GRACE  data  An  independent  test  compared  the  histories  of 
power  inputs  ot£Vs  required  to  explain  measurements  with  those 
of  Pws.  Comparison  results  (Fig.  8)  are  in  excellent  agreement 
Unlike  solar  EUV  radiation  that  most  strongly  heats  the 
thermosphere  at  low-to-mid  latitudes,  stormtime  energy  inputs 
occur  in  the  auroral  oval  and  polar  cap.  Nonetheless,  distinctive 
signatures  of  stormtime  disturbances  appear  at  subauroral 
latitudes  These  include  diminished  TEC  and  |0]/(N2)  ratios  that 
indicate  some  of  the  stormtime  energy  input  is  transported  away 
from  high-latitude  sources.  GRACE  measures  of  mass  density 
profiles  along  the  midnight  meridian  (Fig.  6)  illustrate  the 
equatorward  propagation  of  thermospheric  structures  that  reach 
the  equator  ~~3h  after  energy  injec  tion  at  high  latitudes. 

The  four  main  conclusions  of  our  analysis  of  the  November  10, 
2004  magnetic  storm  are: 

1.  Although  Poynting  flux  distributions  from  DMSP  and  W5 
predictions  sometimes  agree  (Fig.  4),  this  is  not  always  the 
case.  Breakdown  in  agreement  is  most  evident  in  the  midnight 
sector  of  the  auroral  oval.  Integrated  values  of  Poynting  flux 
measured  by  DMSP  agree  fairly  well  with  W5  predictions. 

2.  Energy  added  to  the  stormtime  thermosphere  is  -3  £RC.  A  large 
fraction  of  the  peak  energy  of  ~4  x  10,GJ  reached  the  thermo¬ 
sphere  at  magnetic  latitudes  that  are  not  conjugate  to  the  ring 
current. 

3.  Electromagnetic  power  into  the  thermosphere  during  the 
November  2004  storm,  estimated  as  3£V s  and  Pws.  attained 
peak  values  of  2-3  TW  and  exceed  the  total  EUV  power  to  the 
entire  dayside  of  the  Earth  (Knipp  et  al.,  2005). 

4.  In  agreement  with  dayside  results  of  Bruinsma  and  Forbes 
(2007),  it  takes  about  3  h  after  an  energy  injection  for  a  TAD  to 
reach  the  equator  on  the  nightside. 
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Appendix  A.  :  Volland-Stem  (V-S)  electric  fields 

As  formulated  by  Ejiri  (1978)  the  Volland-Stern  model  offers  a 
simple,  albeit  limited,  method  for  estimating  the  trajectories  of 
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ring-current  ions  in  the  inner  magnetosphere.  One  difficulty  in 
applying  the  model  relates  to  the  connection  between  the 
interplanetary  and  magnetospheric  electric  fields.  Burke  (2007) 
reformulated  V-S  by  combining  it  with  the  Siscoe-Hill  (S-H) 
model  of  the  polar  cap  potential  (0pc)  to  show  that  in  the  absence 
of  shielding  the  electric  field  in  the  inner  magnetosphere  is 
fvs**Pc/2RELv.  where  the  denominator  represents  the  width  of 
the  magnetosphere  along  the  dawn-dusk  line.  This  appendix 
outlines  our  method  for  calculating  eVs  values  used  in  the  main 
text. 

Empirical  studies  suggest  that  on  the  dayside  the  equatorial 
magnetopause  is  nearly  self-similar  in  shape  with  LY ft  1.5  Lx 
(Roelof  and  Sibeck,  1993).  Force  balance  at  the  subsolar  magne¬ 
topause  requires  that 

(M) 

More  precise  estimates  of  magnetopause  dimensions  that 
include  erosion  effects  can  be  introduced  where  more  rigorous 
calculations  are  required. 

Between  1980  and  2000  many  investigators  sought  linear 
relationships  between  0^  and  interplanetary  parameters  (Wei- 
mer,  1995,  Boyle  etal.,  1997).  Burke  et  al.  (1999)  surveyed  analyses 
of  </>pt  dependence  on  solar-wind  parameters  and  showed  that  in 
the  linear  regime  consistently  high  correlation  coefficients  were 
obtained  with 


^E(kV)  =  <*>„  +  f-cVswBr  sin2  |  (A2) 

The  subscript  E  refers  to  interplanetary  electric-field  contribu¬ 
tions;  <P0  is  a  residual  potential  with  typical  values  between  20 
and  30  kV,  £r  =  v  (By2+B/),  and  0  the  IMF  clock-angle  in  the  GSM 
V-Z  plane.  Because  the  clock-angle  enters  IEF  calculations  through 
the  positive-definite  quantity  sin2  0/2,  we  ignore  quadrant 
distributions  and  approximate  6  as  cos  If  we  express 

VswBt  in  k V/Re  (where  1  mV/m  ^ 6.4  kV/RE)  then  the  regression 
slope  Lc  represents  a  3-4  RE  wide  ‘’gate"  in  the  solar  wind  through 
which  geoeffective  streamlines  (equipotentials)  must  pass  to 
reach  the  dayside  magnetopause. 

Siscoe  et  al.  (2002)  suggested  that  during  large  magnetic 
storms  Region  1  currents  generate  perturbation  magnetic  fields 
that  alter  the  shape  of  the  dayside  magnetopause,  thereby 
limiting  access  of  solar-wind  streamlines  and  forcing  the  merging 
rate  to  saturate.  Following  the  suggestion  of  Hill  (1984),  Siscoe  et 
al.  (2002)  postulated  that  a  saturation  potential  <PS  contributes  to 

f*>lK 


(PpC  — 


<pE(ps 

0T+  </>s 


(A3) 


Note  that  when  <PE<z<Ps.  and  the  standard  linear 

relation  Eq  (A2)  is  retrieved.  Siscoe  et  al.  (2002)  derived  an 
expression  for  <PS  in  kV 


where  Psw  is  the  dynamic  pressure  of  the  solar  wind  and  £P  the 
effective  Pedersen  conductance  of  the  polar  ionosphere.  Ober  et 
al.  (2003)  used  data  from  the  SWEPAM  and  MFI  on  ACE  to 
compare  <PFC  predictions  of  S-H  and  other  models  with  DMSP 
measurements  during  the  magnetic  storm  of  March  2001. 
Under  pre-storm  conditions  all  of  the  models  predicted  0PC 
values  measured  by  DMSP.  During  the  main  phase  S-H 
predictions  of  0PC  provided  an  excellent,  upper-bound  envelope 
for  DMSP  observations.  Linear  models  predicted  much  larger 
than  observed  <PPC.  Burke  et  al.  (2007b)  found  that  setting 
=  10  mho  provided  excellent  agreement  between  calculated 


and  observed  0pc  during  a  large  number  of  magnetic  storms 
near  both  equinoxes  and  solstices. 

In  computations  of  £Vs  we  used  hourly  averaged  values  of  the 
solar-wind  density  and  velocity  from  SWEPAM  and  IMF  compo¬ 
nents  in  solar-magnetospheric  coordinates  from  the  MFI  on  ACE. 
In  applications  to  S-H  we  set  (p0  =  25  kV,  Lc  =  3.5,  and  £P  =  10 
mho. 
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